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Abstract 

The cross sections and forward-backward asymmetries of hadronic and leptonic events 
produced in e"^e~ colhsions at centre-of-mass energies from 130 to 183 GeV are presented. 
Results for e"'"e~, t^t~, qq, bb and cc production show no significant deviation from 

the Standard Model predictions. This enables constraints to be set upon physics beyond 
the Standard Model such as four-fermion contact interactions, leptoquarks, Z' bosons 
and R-parity violating squarks and sneutrinos. Limits on the energy scale A of e"'"e~ff 
contact interactions are typically in the range from 2 to 10 TeV. Limits on R-parity 
violating sneutrinos reach masses of a few hundred GeV / for large values of their Yukawa 
couplings. 
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1 Introduction 



In the period 1995-97, the centre-of-mass energy of the LEP e^e~ colhder was increased 
in five energy steps from 130 to 183 GeV, so opening a new energy regime for electroweak 
cross section and asymmetry measurements. Such measurements provide a test of the 
Standard Model (SM) and allow one to place limits on possible extensions to it. 

This paper begins by providing in Section |] the definitions of cross section and 
asymmetry used here. A brief description of the ALEPH detector is given in Section § 
and details of the luminosity measurement and data/Monte Carlo samples in Section ^. 

Section ^ describes the measurement of the qq cross section. It also explores heavy 
quark production, providing measurements of -Rb [Re) which are here defined as the ratio 
of the bb (cc) cross section to the total qq cross section. Constraints on qq forward- 
backward asymmetries are obtained using jet charge measurements. In Section ^, cross 
section and asymmetry results are reported for the three lepton species. 

Based on these results. Section |^ gives limits on extensions to the SM involving contact 
interactions, R-parity violating sneutrinos, leptoquarks, Z' bosons and R-parity violating 
squarks. 

2 Definition of Cross Section and Asymmetry 

Cross section results for all fermion species are provided for 



i) the inclusive process, comprising all events with Js'/s > 0.1, so including events 



ii) the exclusive process, comprising all events with Js'/s > 0.9, so excluding radiative 



events, such as those in which a return to the Z resonance occurs. 

Here, the variable s is the square of the centre-of-mass energy. For leptonic final states the 
variable s' is defined as the square of the mass of the outgoing lepton pair. For hadronic 
final states s' is defined as the mass squared of the Z/7* propagator. This latter choice is 
necessary, because as a result of gluon radiation, (which may occur before or after final 
state photon radiation), the mass of the outgoing quark pair is not well defined. 

Interference effects between ISR and final state radiative (FSR) photons affect the 
exclusive cross sections at the level of a few percent and are not accurately described 
by existing Monte Carlo (MC) generators. They are particularly prominent when 
the outgoing fermions make a small angle to the incoming e"'"e~ beams. To reduce 
uncertainties related to this, the exclusive cross section and asymmetry results presented 
here are defined so as to include only the polar angle region |cos^| < 0.95, where 6 
is the polar angle of the outgoing fermion. This is unnecessary for the inclusive cross 
sections, since they are relatively insensitive to radiative photons. The inclusive results 
are therefore defined to include the full angular acceptance. 

When selecting events experimentally, the variable s'^ is used, which provides a good 
approximation to s' when only one ISR photon is present: 




having hard initial state radiation (ISR). 




sin + sin - |sin(6'i + 6*2)1 



X s . 




sin 6*1 + sin 6*2 + |sin(6'i + 6*2) | 



Here 9i and 62 are the angles of the final state fermions f and f measured with respect to 
the direction of the incoming e~ beam or with respect to the direction of a photon seen 
in the apparatus and consistent with ISR. If two or more such photons are found, the 
angles are measured with respect to the sum of their three-momenta. The fermion fiight 
directions are determined in electron and muon pair events simply from the directions of 
the reconstructed tracks, in tau pair events from the jets reconstructed from the visible 
tau decay products, and in hadronic events from the jets formed when forcing the event 
into two jets after removing isolated, high energy photons detected in the apparatus. 

If an event contains two or more ISR photons, then unless these photons all go in the 
same direction, the variable s'^ ceases to be a good approximation to s' . Such events can 



pass the exclusive selection by being reconstructed with J s'^/ s > 0.9, but nonetheless 



have ys'/s < 0.9. They are called 'radiative background'. 

For dilepton events, the differential cross section is measured as a function of the angle 
6*, defined by 

smi(^+ + ^_) ^ ^ 

where 6- and ^+ are the angles of the negatively and positively charged leptons, 
respectively, with respect to the incoming e~ beam. The angle 6* corresponds to the 
scattering angle between the incoming e~ and the outgoing 1~, measured in the 1"'"1~ rest 
frame, provided that no large-angle ISR photons are present. 

The forward-backward asymmetries are determined from the formula 

^FB = ; 5 (3) 

where ap and as are the cross sections to produce events with the negatively charged 
lepton in the forward < 90°) and backward {6_ > 90°) hemispheres, respectively, 
defined in the same limited angular acceptance as given above. Determining the forward- 
backward asymmetries from this formula avoids any specific assumption about the angular 
dependence of the differential cross section. 



3 The ALEPH Detector 

The ALEPH detector and performance are fully described in |l| and p[. In October 
1995, the silicon vertex detector (VDET) described in these papers was replaced by an 
improved detector P], which is used for the analyses presented here. A brief description 
of the ALEPH detector follows. 

The main tracking detector is a time projection chamber (TPC) lying between radii of 
30 and 180 cm from the beam axis. It provides up to 21 three-dimensional coordinates per 
track. Inside the TPC is a small drift chamber (ITC) and within this, the new VDET. 
The latter has two layers of silicon, each providing three-dimensional coordinates. All 
three tracking detectors contribute coordinates to tracks for polar angles to the beam 
axis up to |cos6'| < 0.95. They are immersed in a 1.5 T axial magnetic field provided 
by a superconducting solenoid. This allows the momentum p of charged tracks to be 
measured with a resolution of o'{p)/p = 6 x 10~^pt © 0.005 (where pt is the momentum 
component perpendicular to the beam axis in GeV/c). The three-dimensional impact 



parameter resolution is measured with an accuracy of (34 + 70/p) x (1 + 1.6cos^6') fim 
(where p is measured in GeV/c). 

Between the tracking detectors and the solenoid is an electromagnetic calorimeter 
(ECAL), which provides identification of electrons and photons, and measures their 
energies E with a resolution of a{E)/E = 0.18/ ^E{GeV) + 0.009. Outside the solenoid, 
is a hadron calorimeter (HCAL), which, combined with the ECAL, measures the energy 

of hadrons with a resolution of a{E)/E = 0.85/ yjE{GeV). The HCAL is also used 
for muon identification, together with muon chambers lying outside it. The ECAL and 
HCAL acceptances extend down to polar angles of 190 and 100 mrad to the beam axis, 
respectively. 

The luminosity is measured with a lead/proportional-chamber electromagnetic 
calorimeter (LCAL) covering the small angle region between 46 and 122 mrad from the 
beam axis. A tungsten/silicon electromagnetic calorimeter (SICAL) covering the angular 
range from 24 to 58 mrad is used to provide a cross-check. 

In general, charged tracks are considered good for the analyses presented here if they 
originate within a cylinder of radius 2 cm and length 10 cm, centered at the interaction 
point, and whose axis is parallel to the beam axis. They must also have at least four TPC 
hits, a momentum larger than 0.1 GeV/c and a polar angle to the beam axis satisfying 
Icos^l < 0.95. 

By relating charged tracks to energy deposits found in the calorimeters and using 
photon, electron and muon identification information, a list of charged and neutral energy 
flow particles is created for each event and used in the following analyses. 

4 Data and Monte Carlo Samples and the 
Luminosity Measurement 

The data used were taken at five centre-of-mass energies which are given in Table [^. This 
table also shows the integrated luminosity recorded at each energy point, together with 
its statistical and systematic uncertainties. For analyses relying heavily on the VDET, 

Table 1: Centre-of-mass energies and integrated luminosities of the high energy 
data samples. The two uncertainties quoted on each integrated luminosity 
correspond to its statistical and systematic uncertainty respectively. 



Energy (GeV) 


Luminosity (pb ^) 


130.2 


6.03 ±0.03 ±0.05 


136.2 


6.10 ±0.03 ±0.05 


161.3 


11.08 ±0.04 ±0.07 


172.1 


10.65 ±0.05 ±0.06 


182.7 


56.78 ±0.11 ±0.29 



such as the study of bb production, data collected at 130 and 136 GeV in 1995 are 
discarded, since the new VDET was not fully installed. However, data taken at these two 
energy points in 1997, corresponding to integrated luminosities of 3.30 and 3.51 pb~^, 
respectively, are used. 



The luminosity is measured using the LCAL calorimeter following the analysis 
procedure described in Ref. 0], with a slightly reduced acceptance due to the shadowing 
of the LCAL detector by the SICAL below 59 mrad. The systematic uncertainty on 
the luminosity, which is assumed to be fully correlated between the different energy 
points, includes a theoretical uncertainty, estimated to be 0.25% from the BHLUMI 
Monte Carlo generator. The luminosity measurement is checked by comparison with the 
luminosity measured independently using the SICAL detector. The two are consistent 
within the estimated uncertainties. 

Samples of Monte Carlo events were produced as follows. The generator 
BHWIDE vl.Ol is used for the electron pair channel and KORALZ v4.2 [|] for the 
muon and tau pair channels. Simulation of diquark events relies on production of the 
initial qq system and accompanying ISR photons with either KORALZ or PYTHIA v5.7 
1^. The simulation of FSR and fragmentation is then carried out with the program 
JETSET v7.4 p[. The PYTHIA generator is also used for four-fermion processes such 



as the Z pair and Ze"'"e~ channels. The programs PHOT02 ||TOl, HERWIG v5.9 |[TT| and 
PYTHIA are used to generate the two-photon events. Finally, backgrounds from W pair 
production are studied using the generators KORALW vl.21 and EXCALIBUR |1T3[|. 



5 Hadronic Final States 



Section |5T^ describes the measurement of the qq cross section. Sections and |5]^ study 
heavy quark production, providing measurements of -Rb (Re) which are here defined as 
the ratio of the bb (cc) cross section to the total qq cross section. The measured values of 
Rh (Re) are statistically independent of the qq cross section measurement. Furthermore, 
by measuring these ratios, rather than the bb and cc cross sections, one benefits from 
the cancellation of some systematic uncertainties. Section places constraints upon 
qq forward-backward asymmetries using a jet charge technique applied to b-enriched and 
b-depleted event samples. 



5.1 The Hadronic Cross Section 

The hadronic event selection begins by requiring events to have at least seven good 
charged tracks. The energy flow particles are then clustered into jets using the JADE 
algorithm with a clustering parameter t/cut of 0.008 . Thin, low multiplicity jets with 



an electromagnetic energy content of at least 90% and an energy of more than 10 GeV 
are considered to be ISR photon candidates. The visible mass M^i^ of the event is then 
measured using charged and neutral energy flow particles, but excluding these photons 
and energy flow particles which make an angle of less than 2° to the beam axis. The 
distribution of Mvis is shown in Fig. |1| for data taken at 183 GeV. It is required to be 
more than 50 GeV/c^. 



The inclusive selection makes the additional requirement that J s'^/ s > 0.1. There is 



actually negligible acceptance for events with 0.1 < y s' / s < 0.3 as a result of the cut on 

Mvis. The inclusive cross sections are therefore extrapolated down to \J s'/ s = 0.1 using 
the KORALZ generator. 



The exclusive selection requires \J s'^/ s > 0.9. Here s'^ is determined from the 
reconstructed jet directions, when the event is clustered into two jets, after first removing 
reconstructed ISR photons. For the exclusive selection, these jets are required to have 
|cos^| < 0.95 . The ^Js'^/s distributions for centre-of-mass energies of 130, 161, 172 and 
183 GeV are displayed in Fig. |, together with the expected background. 

For the exclusive process, two additional cuts are then applied. Firstly, Mvis is 
required to exceed 70% of the centre-of-mass energy. This suppresses residual events 
with a radiative return to the Z. Such events can have \J s'^/ s > 0.9 if they emit two or 

more ISR photons. Figure ^ shows Mvis for events with \J s'^/ s > 0.9. The contribution 
from doubly radiative events is clearly seen at low masses. Secondly, when above the 
W+W~ threshold, (i.e. ^/s > 161 GeV), about 80% of W+W~ background is eliminated 
by requiring that the thrust of the event exceeds 0.85 . The thrust distribution is shown 
in Fig. ||. 




Figure 1: Visible mass distribution at a centre-of-mass energy of 183 GeV 
for events having at least seven tracks. The white histogram shows the expected 
signal, whilst the hashed histograms give the expected contributions of the 77 and 
W pair backgrounds. The arrow indicates the cut used in the inclusive selection. 



The selection efficiencies are estimated using the KORALZ Monte Carlo samples. The 
efficiencies are given for all centre-of-mass energies in Table |^. 

The uncertainties in the efficiencies arise from the statistical uncertainties on the Monte 
Carlo and also from a number of systematic effects, which are assessed as described below. 

i) Uncertainties related to the simulation of ISR are estimated from the difference 
in the efficiencies determined from KORALZ and PYTHIA. Only the generator 
KORALZ simulates ISR using YFS exponentiation [§]. 



130GeV ALEPH 161 GeV ALEPH 




0.4 0.6 0.8 1 0.4 0.6 0.8 1 



y/s-J s Vs'„/s 

Figure 2: \/ s'^/ s distribution for hadronic events at centre-of-mass energies from 
130 to 183 GeV. This is compared to the Monte Carlo expectations shown by the 
white histograms. The hashed areas correspond to background contributions and 
are dominated by W pair production. The distribution at 136 GeV is omitted 
because it closely resembles that at 130 GeV. 
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Figure 3: Visible mass distribution at 183 GeV centre-of-mass energy for 
hadronic events with \/s'^/s > 0.9. This is compared with Monte Carlo 
expectations shown by the white histogram. The hashed area represents the 
expected contribution from events with generated -\/ s'/s < 0.9. The arrow 
indicates the cut used in the exclusive selection. 
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Figure 4: Thrust distribution at 183 GeV centre-of-mass energy for hadronic 
events with ^s!^/s > 0.9. The white histogram shows the expected signal and 
the hashed one corresponds to the W pair background. The arrow indicates the 
cut used in the exclusive selection. 



Table 2: Selection efficiencies and background fractions. For e"'"e production, 
numbers are given for (1) -0.9 < cos 6'* < 0.9 and (2) -0.9 < cos 6** < 0.7. 



Js'/s 
V ' 




Event 


Efficiency 


Background 


cut 


(GeV) 


type 


(%) 


(%) 


n 1 


1 SO 

-LOW 


nri 
HH 


SQ 6 + Q 


2 4 ± 3 








1 O . O _l_ -L . \J 


4 ± 1 








55.3 ± 1.1 


5.7 ± 2.4 






nn 
HH 


on c _|_ n 7 


1 « _l_ f) 

X 1 yj .Zj 






1 1 1 

LI LI 


70 O _|_ 1 c 


^ ± 1 








53.9 ± 1.3 


8.9 ±4.0 




161 




88.4 ±0.6 


4.0 ±0.1 






r r 


76.5 ± 1.5 


3.3 ± 1.2 






T+T- 


11.2 ± 1.1 


8.3 ± 1.0 




172 


qq 


87.3 ±0.6 


10.2 ±0.1 








75.6 ± 1.6 


4.1 ± 1.3 








44.8 ± 1.0 


11.3 ±3.7 




183 


qq 


83.9 ±0.7 


17.8 ±0.2 








75.6 ± 1.6 


5.4 ± 1.2 






f+f- 


47.4 ± 0.9 


12.9 ±2.6 


0.9 


130 


HH 


92.2 ± 1.0 


9.5 ± 0.3 
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136 


no 


89 2 ± 8 


8.8 ± 0.4 
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161 


no 

HH 


89 3 ± 7 


7.7 ± 0.3 






e e i^ij 


on 1 _(_ 1 p: 


1 n n _(_ n q 








q4 -L 1 Q 


11 n + 1 q 

X X . W _l_ X . 






1 1^ 1 1 

LI LI 


Qfi + 1 8 

C/U.U X .0 


4 c _l_ n q 






T~^T 
/ / 


«0 Q _|_ 1 7 

W.^ . (J _l_ X.I 


q _c: _|_ 7 




172 




90.2 ± 0.6 


6.4 ± 0.3 






e+e-(l) 


90.7 ± 1.5 


10.7 ± 1.0 






e+e-(2) 


96.3 ± 1.2 


11.9 ± 1.4 








96.8 ± 1.8 


6.0 ± 1.3 






T+T- 


61.6± 1.7 


13.9 ± 3.8 




183 


qq 


85.9 ±0.8 


9.2 ±0.1 






e+e-(l) 


87.5 ± 1.4 


10.9 ±0.7 






e+e-(2) 


95.2 ± 1.2 


12.4 ±0.9 






//+//- 


95.9 ± 1.8 


6.5 ± 1.3 






-+-- 


67.9 ± 1.8 


13.3 ± 3.0 



ii) The overall ECAL energy calibration is determined each year, from Bhabha events 
recorded when running at the Z peak. Its uncertainty is estimated to be ±0.9%, 
based upon a comparison of the detector response to low energy 1 GeV) 
electrons in 77 —>■ e~^e~ events, with that to high energy electrons in Bhabha events. 
The calibration of the HCAL energy scale uses minimum ionizing particles. The 
statistical uncertainty on this calibration is ±2%. These uncertainties on the 
calorimeter energy scales are considered to be uncorrelated from year to year. 

iii) The energy response to 45 GeV jets, in data and Monte Carlo events produced at the 
Z resonance, is compared and the difference parametrized as a function of polar angle 
to the beam axis. The selection efficiency is corrected using this parametrization. 
The systematic uncertainty due to the energy response is derived from the change in 
the measured cross section when this correction is applied. This uncertainty, which 
is the largest of the detector related ones, is considered to be correlated from year 
to year. 

iv) The measured polar angles of jets with respect to the beam axis can suffer from 
small systematic biases, particularly at low polar angles. These biases are studied 
using events produced at the Z resonance. Differences between the data and Monte 
Carlo are taken into account via their effect on s'^. The associated systematic 
uncertainty is given by the change in the measured cross section when applying 
these corrections. It is considered to be correlated from year to year. 

The background fractions remaining after the selection are also given in Table ^. The 
uncertainty in these backgrounds receives contributions from the detector calibration and 
energy response, estimated as above. Additional uncertainties are studied as follows: 

i) The two-photon background is simulated with the PYTHIA, PHOT02 and HERWIG 
generators. It is normalized to the data in the region Mvis < 50 GeV/c^ (Fig. 0), and 
the difference between the expected cross section from PIIOT02 and the normalized 
one is used to calculate the uncertainty on the process. This background dominates 
at centre-of-mass energies of 130 and 136 GeV. 

ii) Above the W pair production threshold, the W'*'W^ background is estimated 
using KORALW. At 161 GeV the input cross section is taken from the theoretical 



expectation, computed with the GENTLE |jT5[ program for a W mass of 80.39 ± 



0.06 GeV/c^|T^. At 172 and 183 GeV the cross sections are taken from the average 
measurement of the four LEP experiments |T^. Uncertainties in the W pair cross 
sections are propagated to estimate the associated systematic uncertainties on the 
qq cross sections. 

iii) Other four-fermion backgrounds such as ZZ and Ze"'"e~ are estimated using PYTHIA 
and EXCALIBUR. They introduce only a ±0.1% systematic uncertainty on the 
exclusive qq cross section measurement. 

iv) Events where ISR occurs from the incoming electron and positron remain an 
important background for the exclusive selection. Systematic uncertainties on this 
background arise from the simulation of ISR photons and from the energy scale and 
response of the detector. These are already taken into account as described above. 



The cross section measurements at each centre-of-mass energy are hsted in Table 
where they are compared with the SM predictions. These predictions are discussed in 
Section [7.1|. For the exclusive results, the measurements and predictions refer to qq final 



states with |cos6'| < 0.95, where d is the polar angle of the quark with respect to the beam 
axis. A breakdown of contributions to the systematic uncertainties on the cross section 
measurements is given in Table ^. 

5.2 Measurement of the bb Production Fraction i?b 

The ratio -Rb of the bb to qq production cross sections is measured at centre-of-mass 
energies in the range 130-183 GeV. Hadronic events are chosen following the exclusive 



selection described in Section ^T^- The bb events are separated from the hadronic ones 
using the relatively long lifetime of b hadrons. From the measured impact parameters of 
the tracks in an event, the confidence level that all these tracks originate from the primary 



vertex is calculated and required to be less than a certain cut, which is chosen to 
minimize the total uncertainty on (7^5 • 

The estimated selection efficiency and background fraction at each centre-of-mass 
energy are given in Table ^ At low centre-of-mass energies about two thirds of the 
background is due to radiative events, with the remainder being dominated by cc events. 
At higher energies, radiative background, cc events and background from four-fermion 
events contribute in roughly equal proportions. 

The largest contribution to the systematic uncertainty arises from the b tagging 
efficiency. This is evaluated by using the same technique as above to measure the fraction 
i?b of bb events in Z data taken in the same year. The difference between this measurement 
and the very precisely known world average measurement of i?b at the Z peak is then taken 
as the systematic uncertainty. The b tagging efficiency is almost independent of the centre- 
of-mass energy, because the impact parameters of tracks from a decaying b hadron depend 
little on its energy. The systematic uncertainties evaluated at the Z peak can therefore 
be directly translated to higher centre-of-mass energies. The small decrease in tagging 
efficiency with centre-of-mass energy seen in Table ^ is due to a decay length cut, which 
the b tag uses to reject tracks from K° decay. This introduces no significant additional 
uncertainty on the tagging efficiency. 

Systematic uncertainties due to the charm and light quark backgrounds are smaller. If 
one neglects the dependence of the quark production fractions on centre-of-mass energy, 
they are already taken into account by the study using Z peak data mentioned above. The 
impact parameter resolution is monitored using tracks with negative impact parameters. 

The number of observed events, together with the measured and predicted values of 
i?b are given in Table Figure |^ summarizes the measurements of i?b as a function of 

5.3 Measurement of the cc Production Fraction Rc 

The ratio -Rc of the cc to qq production cross sections is measured at a centre-of-mass 
energy of 183 GeV. Hadronic events are chosen using the exclusive selection described 



in Section |5.1| , and clustered into two jets using the JADE algorithm. An additional 
acceptance cut requiring that both jets have |cos^| < 0.9 is applied to ensure that the 
event is well contained inside the VDET. 



Table 3: Measured cross sections with statistical and systematic uncertainties 
for different channels at centre-of-mass energies from 130 to 183 GeV. The SM 
predictions are also given, together with the number of selected events (before 
background subtraction). The exclusive cross sections correspond to the restricted 
angular range |cos^| < 0.95. For the Bhabha process, results are given for (1) 
-0.9 < cos 6** < 0.9 and (2) -0.9 < cos 6** < 0.7. 



/ // 


-t-'cm 


bvent 


JNo. 


(7 


O TV yr „ J • „ J. • „ 

bM prediction 


cut 


(GeV) 


type 


Events 


(Pb) 


(pb) 


0.1 


130 


qq 


1858 


335.6 ±7.9 ±4.5 


327.0 








110 


22.5 ±2.2 ±0.5 


21.9 








94 


25.9 ±2.9 ±0.6 


21.9 




136 


qq 


1558 


280.8 ±7.2 ±3.7 


269.8 








103 


20.4 ±2.1 ±0.5 


18.7 








67 


17.8 ±2.5 ±0.5 


18.6 




161 


qq 


1520 


149.0 ±4.0 ±1.7 


147.6 








107 


12.2 ± 1.2 ± 0.3 


11.2 








78 


14.6 ±1.8 ±0.3 


11.2 




172 


qq 


1270 


122.6 ± 3.8 ± 1.2 


122.8 








74 


8.8 ±1.1 ±0.2 


9.5 








51 


9.5 ± 1.5 ±0.3 


9.5 




183 


qq 


6072 


104.8 ±1.6 ±0.9 


104.4 








406 


8.84 ±0.47 ±0.19 


8.22 






1 


2 11 


7.80 ± 0.58 ± 0.17 


8.22 


0.9 


130 


qq 


440 


71.6 ±3.8 ±1.1 


70.7 






e+e-(l) 


1186 


191.3 ±6.2 ±3.5 


186.7 






e+e-(2) 


274 


41.1 ±2.8 ±0.9 


39.0 








48 


7.9 ±1.2 ±0.2 


7.0 








49 


10.9 ±1.8 ±0.4 


7.3 




136 


qq 


351 


58.8 ±3.5 ±0.9 


57.3 






e+e-(l) 


1051 


162.2 ±5.6 ±3.5 


167.3 






e+e-(2) 


212 


29.5 ±2.4 ±0.7 


33.9 








43 


6.9 ± 1.1 ±0.2 


6.1 








27 


5.6 ±1.3 ±0.2 


6.3 




161 


qq 


321 


29.94 ±1.8 ±0.4 


30.7 






e+e-(l) 


1393 


119.7 ±3.7 ±2.3 


119.0 






c+c-(2) 


302 


25.6 ± 1.7 ±0.5 


24.8 








50 


4.49 ± 0.69 ± 0.09 


3.88 








44 


5.75 ±0.96 ±0.23 


4.01 




172 


qq 


271 


26.4 ± 1.7 ± 0.4 


25.1 






c+c-(l) 


1166 


107.8 ±3.5 ±2.1 


102.5 






e+e-(2) 


268 


23.0 ± 1.6 ±0.5 


21.5 








29 


2.64 ±0.53 ±0.06 


3.32 








26 


3.26 ±0.74 ±0.09 


3.43 




183 


qq 


1165 


21.71 ±0.70 ±0.23 


21.1 






c+c-(l) 


5063 


90.9 ±1.4 ±1.7 


90.9 






e+e-(2) 


1171 


18.99 ±0.63 ±0.36 


19.11 








175 


2.98 ±0.24 ±0.06 


2.89 






^+^- 


129 


2.90 ±0.29 ±0.09 


2.98 



Table 4: Contributions to the systematic uncertainties on the qq cross section 
measurements, for all energies and for both inclusive and exclusive processes. All 
quoted values are in percent. 





Description 




E 


cm (GeV) 




cut 




130 


136 


161 


172 


183 


0.1 


MC statistics 


0.3 


0.3 


0.3 


0.3 


0.1 




ISR simulation 


0.3 


0.3 


0.3 


0.3 


0.3 




Energy scale 


0.4 


0.5 


0.4 


0.4 


0.3 




Detector response 


0.6 


0.5 


0.5 


0.5 


0.6 




77 qq 


0.3 


0.2 


0.05 


0.05 


0.04 




w+w- 






0.05 


0.05 


0.04 




Luminosity 


1.0 


1.0 


0.7 


0.7 


0.5 


0.9 


MC statistics 


0.4 


0.4 


0.4 


0.4 


0.2 




ISR simulation 


0.3 


0.4 


0.4 


0.7 


0.4 




Energy scale 


0.4 


0.3 


0.3 


0.3 


0.3 




Detector response 


0.9 


0.9 


0.7 


0.7 


0.7 




W+W" 






0.02 


0.02 


0.01 




ZZ 






0.01 


0.01 


0.01 




Other four-fermion 






0.03 


0.03 


0.03 




Luminosity 


1.0 


1.0 


0.7 


0.7 


0.5 



Table 5: Selection efficiencies and background fractions for the measurement of 



E 

-'-'cm 


Efficiency 


Background 


(GeV) 


(%) 


(%) 


130 


40.3 ± 1.0 


12.1 ± 1.6 


136 


40.8 ± 1.0 


10.0 ± 1.5 


161 


40.7 ± 1.5 


10.8 ± 1.2 


172 


38.0 ± 1.4 


10.6 ± 1.3 


183 


37.8 ±0.9 


9.6 ± 1.0 



Table 6: Measured values of i?b with statistical and systematic uncertainties for 
bb production with \/ s' / s > 0.9 and \cos9\ < 0.95, where 9 is the polar angle of 
the b quark with respect to the beam axis. The SM predictions are also given, 
together with the number of selected events (before background subtraction). 



E 


No. 




SM prediction 


(GoV) 


E\xaits 






130 


19 


0.176 ±0.044 ±0.004 


0.190 


136 


20 


0.214 ±0.050 ±0.005 


0.186 


161 


24 


0.159 ±0.034 ±0.006 


0.175 


172 


16 


0.134 ±0.036 ±0.007 


0.173 


183 


91 


0.176 ±0.019 ±0.005 


0.171 
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Figure 5: Measured values of the ratio i?b at various centre-of-mass energies, 
compared with the SM expectations. 



Background from bb events is suppressed by taking advantage of the relatively long 
lifetime and high mass of b hadrons. Three algorithms are used sequentially. The first 
rejects events on the basis of track impact parameter significances the second using 



the decay length significance of reconstructed secondary vertices [0 and the third relying 



on a comparison of the total invariant mass of high impact parameter significance tracks 



with the charm hadron mass IT^. Together the hadronic selection and bb rejection have 



a 60% efficiency for cc events with ys'/s > 0.9. The residual fraction of bb events in the 
resulting sample is 4%. 

The final cc selection uses a neural network. This was trained to separate cc jets 
(giving a neural network output close to one) from light quark jets (giving an output 
close to minus one). This network uses twelve variables per jet. These are listed below, 
ordered according to decreasing weight. 

• The sum of the rapidities with respect to the jet axis of energy fiow particles within 
40° of this axis. 

• The sphericity of the four most energetic energy fiow particles in the jet, calculated 
in their rest frame. 

• The total energy of the four most energetic energy fiow objects in the jet. 

• The number of identified leptons (electrons or muons) in the jet with momentum 
larger than 1.5 GeV/c. 

• The transverse momentum squared p']_ with respect to the jet axis of the TTgoft 
candidate from D* — > vTsoftX, defined as the charged track in the jet with the smallest 
value of p± and an energy between 1 and 4 GeV. 

• The confidence level that all charged tracks in the jet originate from the primary 
vertex. 

• The energy of a sub jet of mass 2.1 GeV/c^ built around the leading energy fiow 
particle in the jet. 

• The momentum of the leading energy fiow particle in the jet. 

• The number of energy fiow objects in a cone of half-angle 40° around the jet axis. 

• The confidence level that all charged tracks in the jet having a rapidity with respect 
to the jet axis exceeding 4.9, originate from the primary vertex. 

• The decay length significance of a reconstructed secondary vertex. 

• The energy of a reconstructed D meson (if any). D meson reconstruction is 
attempted in the charged decay channels Kvr, Kirn and KTrvrvr. 

The distribution of the sum of the neural network outputs for the two jets in each event 
is shown in Fig. ^. A lower cut is placed on this to reject light quark events. It was also 
found that an upper cut served to reject a tail of remaining bb events. The final selection 
efficiencies and event fractions for each fiavour are summarized in Table ^. 

To evaluate the systematic uncertainty introduced by a given selection variable, the 
ratio of distributions of this variable in data to Monte Carlo is determined. After 
smoothing to reduce statistical fiuctuations, this ratio is then used to reweight individually 
bb, cc and light quark event fiavours in the Monte Carlo. The relative change in the 
resulting efficiency for each fiavour is taken as its systematic uncertainty. A cross-check 
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Figure 6: Distribution of the sum of the neural network outputs for hadronic 
events at 183 GeV centre-of-mass energy after hadronic selection and bb rejection. 
The lines indicate the chosen cuts. 



is performed by measuring Rc on the Z data taken in the same year; the measured value 
is consistent with the SM prediction. All contributions to the systematic uncertainty are 
given in Table ^. 

In the data, 153 events are selected, of which 74.8 are estimated from the Monte Carlo 
simulation to be background. The value of Rc at ^/s = 183 GeV, for ^ s' / s > 0.9 and 
|cos0| < 0.95, where 6 is the polar angle of the c quark with respect to the beam axis, is 
found to be 

Rc = 0.276 ±0.041(stat) l°:°i^(syst) - 0.284 x (Kds - 0.584) - 0.390 x (i?b - 0.171) , (4) 

where the dependence on a possible departure of the light quark or bb fractions from their 
SM expectations is explicitly given. The measured value of Rc is consistent with its SM 
expectation of 0.244 . 

5.4 Measurement of A^g using Jet Charge 

Constraints upon the forward-backward asymmetries Ap^ of qq events at y/s = 183 GeV 
are obtained using a jet charge technique. This is performed separately for b-enriched 
and b-depleted events, which allows the asymmetry in bb events to be well determined. 



Table 7: Selection efficiency 

for each flavour, and its expected 
fraction in the final event sample, 
assuming SM cross sections. 





Efficiency (%) 


Fraction (%) 


uu 


6.18 ±0.25 


14.8 ±0.6 


dd 


6.66 ±0.33 


10.3 ±0.5 


ss 


6.09 ±0.31 


9.7 ±0.5 


cc 


22.07 ±0.48 


51.8 ±0.8 


bb 


8.61 ±0.37 


13.4 ±0.6 



Table 8: Contributions to the 

systematic uncertainty on the measured 
value of Rn. 



Description 


Uncertainty 


Luminosity 


±0.001 


Efficiency of qq selection 


±0.001 


Background in qq selection 


±0.001 


Performance of bb rejection tag 


+0.012 
-0.003 


Performance of neural network 


+0.006 
-0.010 


MC statistics 


±0.007 


Total 


+0.016 
-0.013 



The analysis uses the hadronic events with ys'/s > 0.9, selected as described in 
Section ^TT|. To ensure that the events are well contained in the tracking chambers, they 
are also required to satisfy |cos0*| < 0.9. The events are then divided into two samples, 
according to whether they pass or fail the b lifetime tag of Section (here used with 
looser cuts). This gives a 91% pure sample of bb events plus a sample of predominantly 
light / charm quark events. 

After clustering all the events into two jets, the jet charge Qjet of each jet is determined, 
where 

■^track / -^track 

Qiet= E P\\^Q^/ E (5) 
i=l / i=l 

and the sums extend over the charged tracks in the jet. The track momentum component 
parallel to the jet and its charge are p\\i and Qi, respectively. The parameter k is set to 0.3, 
because this minimizes the uncertainty on the final result. The mean charge difference 
between the forward and backward jets {Qfb) = {Qjct) ~ (Qfct) is then formed. The ffist 
row of Table ^ shows the observed value of (Qfb) in the data. 



Table 9: Comparison of (Qfb) in data with the prediction of Equation |6| using 
SM values of cjq and Ap^ . Also shown is the contribution to this prediction from 
signal and background events. 



{Qfb) 


b tagged events 


anti-b tagged events 


Data 
Expectation in SM 
Contribution from signal 
Contribution from background 


-0.029 ±0.018 
-0.052 
-0.048 
-0.004 


0.022 ±0.007 
0.014 
0.018 
-0.003 



The expected value of (Qfb) depends on the qq cross sections cXq and asymmetries 
Ap-Q. Defining as the selection efficiency for event type i, the prediction is 



(Qfb) = -^^ , (6) 

q X 

where the sums extend over the quark flavours q and the various background types x. 
The parameters 6^ = (Qfet) ~ {Qfet) S^^^ ^^e mean charge separation between the jet 
containing the quark and that containing the antiquark. Table |10| shows the parameters 
6q for each flavour. They are obtained from Monte Carlo simulation, but with additional 
corrections based upon a precise analysis at the Z peak |T^. The parameters 6q vary 
only slowly with centre-of-mass energy so these corrections are still applicable at LEP2. 
The most important correction is that to 6c, which amounts to 26%. This is due to an 
inadequate simulation of charm hadron decay modes. However, even this only alters the 
predicted value of (Qfb) by about one half of the statistical error in the data. Systematic 
uncertainties on the 6q can therefore at present be neglected. The parameters 0.98) 
give the dilution in (QFB)q = ^FB'^q-^q caused by the small angular dependence of the 
acceptance efficiency. 



Table 10: Mean jet charge separations for each quark flavour. 



5n 


5d 


5s 




5b 


0.205 


-0.130 


-0.153 


0.155 


-0.108 



Assuming SM values of cXq and ^fb^ ^^e predicted values of (Qfb) are given in the 
second row of Table ^. They agree with those measured in the data, both for the b- 
enriched and the b-depleted samples. By comparing the measured value of (Qfb) in 
each of these two samples with the predictions of Equation ^, two independent constraint 
equations are obtained for the allowed values of cXq and Ap-g. 

Providing that deviations from the SM are small, it is convenient to approximate each 
of these two constraints by a linear equation. These two linear equations will be used 
when placing limits on physics beyond the SM in Section 0. 

If (Tq and differ from their SM predictions by Aa^ and AAp^, respectively, then 
a Taylor expansion of Equation ^ yields 



/n \ /n \ '9(Qfb) ^ ^ , ^(Qfb) ^ ^q 

{Qfb) - (<;?fb)sm = ^^q + g^q^ ^^fb 



(7) 



Dividing this equation throughout by the estimated uncertainty on the measured value 
of (Qfb) in the data gives an equation of the form 

7 = 5]«qAaq+/5qAA^B . (8) 

q 

Here 7 is the difference between the measured value of (Qfb) and that expected in the 
SM, divided by the measurement uncertainty (such that 7 has an uncertainty of ±1). The 
coefficients and /3q, defined by this equation, can be expressed in terms of e^, 5q, Dq 
and (Qfb)^- The values of 7, and /5q for the b-enriched and b-depleted event samples 



are given in Table 11. The uncertainties in Oq and Pq can be neglected. The equation 



obtained from the b-enriched sample can be interpreted as a measurement of ^fb; since 
Pb dominates the other coefficients. It gives a measured value of A^-q = 0.33 ± 0.19, 
compared with the SM prediction of 0.57 . 



Table 11: Coefficients of the hnear constraint Equation ^ derived for the b- 
enriched and b-depleted data samples. The units of a are pb~^. 



Event Sample 


7 


au 


ad 


as 




at 




Pd 




Pc 


Pb 


b tagged 


1.3 ±1.0 


0.0 


0.0 


0.0 


0.0 


-0.1 


0.0 


0.0 


0.0 


0.4 


-5.4 


anti-6 tagged 


1.2 ± 1.0 


0.7 


-0.6 


-0.7 


0.5 


-0.2 


6.6 


-2.7 


-3.4 


5.1 


-1.1 



6 Leptonic Final States 
6.1 The e+e" Channel 

Electron pair events are selected by requiring the presence of two good tracks of opposite 
charge with a polar angle to the beam axis of |cos ^| < 0.9. The sum of the momenta of the 
two tracks must exceed 30% of the centre-of-mass energy. The total energy associated with 
them in the ECAL must be at least 40% of the centre-of-mass energy. When calculating 
this energy, if one of the tracks passes near a crack in the ECAL, then the associated HCAL 
energy is included if it matches the track extrapolation within 25 mrad. Furthermore, 
the energy of bremsstrahlung photons is included when within cones of 20° around each 
track. 



Only the exclusive cross section is measured, so the requirement ys^/s > 0.9 is 
applied. The main background is due to ISR and this is reduced to a level of 10 - 12%, by 
requiring that the invariant mass of the e^e~ final state, determined from the measured 
track momenta, exceeds 80 GeV/c^. The resulting invariant mass distribution is shown in 
Fig. ^ The discrepancy between data and Monte Carlo seen in this figure arises from an 
inaccurate simulation of the momentum resolution of tracks at small angles to the beam 
axis. As the cut applied on the electron pair invariant mass is very loose, this does not lead 
to a significant systematic uncertainty. The normalization of the radiative background is 
determined from a fit to Fig. |^, using the expected shapes of the signal and background. 
The statistical uncertainty on the fit result leads to the systematic uncertainty due to 
radiative background quoted in Table |l^. 

The selection efficiencies are determined from BHWIDE Monte Carlo samples at each 
energy point. They are given, together with the estimated background, in Table 0. The 
uncertainties on the efficiencies are estimated by varying the calorimeter energy scales as 



in Section 5.1 



The exclusive cross section is determined in two polar angle ranges: —0.9 < cos 9* < 
0.9 and —0.9 < cos 6'* < 0.7. It is dominated by t channel photon exchange, particularly 
in the forward region. The cross section measurements are given in Table ^. The 
contributions to the systematic uncertainties on these measurements are given in Table |T^. 
Uncertainties from possible bias in the polar angle measurement of tracks are negligible. 
They have been assessed by redetermining the cross section without using tracks having 
Icos^l > 0.9. 
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Figure 7: Electron pair invariant mass distribution at 183 GcV ccntrc-of-mass 
energy. The data, which is shown by points, is compared with the Monte Carlo 
expectations for the signal and the radiative background. 



Table 12: Contributions to the systematic uncertainties on the e+e~ exclusive 
cross section measurements, for both cos^* intervals. All quoted values are in 
percent. 





cos^* range 


Description 


Ecm (GeV) 


cut 






130 


136 


161 


172 


183 


0.9 


-0.9 < cosr < 0.9 


MC statistics 


0.4 


0.4 


0.4 


0.4 


0.4 






Energy scale 


1.7 


1.6 


1.4 


1.6 


1.6 








0.01 


0.01 


0.01 


0.01 


0.01 






Radiative background 


0.9 


1.0 


0.9 


1.0 


0.9 






Luiiiinosit}' 


1.0 


1.0 


0.7 


0.7 


0.5 


0.9 


-0.9 < cosr < 0.7 


MC statistics 


0.6 


0.5 


0.6 


0.6 


0.6 






Energy scale 


1.4 


1.3 


1.2 


1.2 


1.1 








0.2 


0.04 


0.03 


0.03 


0.03 






Radiative background 


1.3 


1.6 


1.3 


1.4 


1.3 






Luminosity 


1.0 


1.0 


0.7 


0.7 


0.5 



6.2 The fi+fi- Channel 



The muon pair selection requires events to contain two good, oppositely charged tracks 
with momenta exceeding 6 GeV/c and angles to the beam axis of |cos^| < 0.95. The 
scalar sum of the momenta of the two tracks must exceed 60 GeV/c. The total number of 
good charged tracks in the events must be no more than eight. To limit the background 
from cosmic ray events, both tracks are required to originate near the primary vertex, 
and to have at least four associated ITC hits, which confirms that they were produced 
within a few nanoseconds of the LEP beam crossing. 

Both tracks must be identified as muons, where muon identification is based either on 
the digital hit pattern associated with a track in the HCAL or on its energy deposition 
in the calorimeters: 



The track should fire at least 10 of the 23 drift tube planes in the HCAL. It should 
also fire at least half of the HCAL planes which the track is expected to cross (taking 
into account HCAL cracks), and furthermore should fire 3 or more of the outermost 
10 planes of the HCAL. Alternatively, the track should have at least one associated 
hit in the muon chambers. 

The energy deposition must be consistent with that of a minimum ionizing particle, 
i.e., the sum of the energies associated with the track in the ECAL and HCAL should 
not exceed 60% of the track momentum. Moreover, the sum of these energies and 
the track momentum should be smaller than 60% of the centre-of-mass energy. To 
control a small misidentification background related to calorimeter cracks, tracks 
are also required to have at least one associated hit in the 10 outermost layers of 
the HCAL. 



In the case of the exclusive selection, ys'^/s >0.9, it is also required that the muon 
pair invariant mass exceed 110 GeV/c^, reducing background from radiative events by 
about 40%. The invariant mass distribution is shown in Fig. 

The efficiency of the kinematic selection cuts is estimated using KORALZ Monte Carlo 
events, whilst the muon identification efficiency is measured using muon pair events in 
data recorded at the Z peak in the same year. This efficiency is typically uncertain by 
about ±1.6% as a result of the limited number of Z events. 

For the inclusive process, the main background contamination stems from 77 fi~^fi~ . 
The systematic uncertainty associated with the normalization of this background is 
estimated by comparing data and Monte Carlo in the region of /i'^yU" mass below 
50 GeV/c^, which is shown in Fig. ^ 

For the exclusive process, the main background comes from radiative events and is 
assessed using the mass region below O.Q-y/i of Fig. Residual cosmic ray contamination 
is estimated by relaxing cuts on track impact parameters, and is found to be less than 
0.5%. 

Efficiencies and background levels are summarized in Table ^ and the cross section 
measurements are given in Table ^. The contributions to the systematic uncertainties on 



the cross sections are given in Table 13. 
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Figure 8: Distribution of the invariant mass of f^^fJ-" events at ^/s = 183 GeV, 
with yj s'^l s > 0.9. The data points are compared to the white histogram of 
Monte Carlo \Jir events, normalized to the same integrated luminosity. The 
multi-radiative background is shown as the shaded histogram. The vertical hue 
shows the cut used for the exclusive selection. 



Table 13: Contributions to the systematic uncertainties on the muon pair 
cross section measurements, for all energies and for both inclusive and exclusive 
processes. All quoted values are in percent. 





Description 


Eem (GeV) 


cut 




130 


136 


161 


172 


183 


0.1 


MC statistics 


0.4 


0.4 


0.4 


0.4 


0.4 




\i identification 


2.0 


2.0 


1.9 


2.0 


1.9 






0.0 


0.0 


1.0 


1.0 


0.7 






0.04 


0.05 


0.1 


0.1 


0.1 




w+w- 






0.05 


0.2 


0.3 




Luminosity 


1.0 


1.0 


0.7 


0.7 


0.5 


0.9 


MC statistics 


0.4 


0.4 


0.4 


0.3 


0.4 




11 identification 


1.9 


1.9 


1.8 


1.8 


1.9 






0.00 


0.00 


0.02 


0.1 


0.1 




w+w- 






0.0 


0.03 


0.2 




Radiative background 


0.5 


0.5 


0.6 


0.8 


0.6 




Luminosity 


1.0 


1.0 


0.7 


0.7 


0.5 



40 



> 120 

o 
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Figure 9: Distribution of the invariant mass of fi^ events at 183 GeV centre- 
of-mass energy in the low mass region where the process 77 — > /i"*"/^" is dominant. 
The Monte Carlo prediction is indicated by the shaded histogram. 



6.3 The T+T- Channel 

The tau pair selection begins by clustering events into jets, using the JADE algorithm with 
the clustering parameter ?/cut equal to 0.008 . Tau jet candidates must contain between 
one and eight charged tracks. Events with two tau jet candidates are selected, providing 
that the invariant mass of the two jets exceeds 25 GeV/c^. This requirement removes a 
large part of the 77 background. 



Following an approach already used for tau identification at LEPl ||2^, each tau jet 
candidate is analysed and classified as a tau lepton decay into an electron, a muon, charged 
hadrons or charged hadrons plus one or more tt*^. Both tau jets must be classified in this 
way. Events are required to have at least one tau jet candidate identified as a decay into a 
muon or into charged hadrons or charged hadrons plus 7r°. To suppress background from 
77 — i> ji^ ji~ or events with two jets classified as muonic tau decay are discarded. 

For the exclusive selection, W pair events in which both W decay leptonically represent 
an important background. However, most of this background is rejected by requiring that 
the acoplanarity angle between the two taus be less than 250 mrad. 

Estimated selection efficiencies and background levels are given in Table |^. A more 
detailed breakdown of the selection efficiencies for the various tau decay channels is given 



in Table |14|. The main uncertainty on the selection efficiency arises from the energy 
scale of the calorimeters and is estimated as in Section |5.1| . For the inclusive cross 
section the dominant background comes from 77 t^t^ . The systematic uncertainty 
associated with the normalization of this background is estimated from a comparison 
of data and Monte Carlo in the low visible mass range of selected tau pair events 



(15 GeV/c^ < Mvis < 50 GeV/c^). In the exclusive selection, the dominant background 
is Bhabha events, which sometimes pass the selection criteria if they enter cracks in the 
ECAL acceptance. 

The cross section measurements are listed in Table ^ The contributions to the 



systematic uncertainties are given in Table 15 



Table 14: Percentage efficiency of the t~^t selection for the various tau decay 
channels at a centre-of-mass energy of 183 GeV. 



Ti decay 


/ 


T2 decay 


^Js'/s >0.1 


^s'/s 


>0.9 




/ 




0.8 ± 0.5 


1.1 


± 1.1 




/ 


hadrons Ur 


56.3 ± 1.8 


77.1 


± 2.5 




/ 


hadrons tt^ Vr 


61.7 ± 1.3 


87.3 


± 1.5 




/ 




49.2 ± 2.1 


74.9 


± 2.9 


hadrons 


/ 


hadrons 


48.0 ± 2.2 


67.5 


± 3.3 


hadrons I'r 


/ 


hadrons 7r° 


53.6 ± 1.2 


77.9 


± 1.6 


hadrons z/^ 


/ 




31.6 ± 1.7 


50.2 


± 2.9 


hadrons tt^ Vr 


/ 


hadrons 7r° Ur 


58.7 ± 1.2 


79.9 


± 1.6 


hadrons ti^ 


/ 




35.2 ± 1.3 


48.5 


± 2.1 




/ 




0.0 


0.0 





6.4 Measurement of the Lepton Asymmetries 

Figures 0, [11] and |12| show the observed cos 6* distributions for electron, muon and tau 
pair events passing the exclusive selections at a centre-of-mass energy of 183 GeV. 

As discussed in Section 0, the dilepton asymmetries are determined from the fraction 
of events in which the negatively charged lepton enters the forward/backward hemispheres 
with respect to the incoming electron. They are measured only for the exclusive process 
and defined in the range |cos^| < 0.95, where 6 is the polar angle of the outgoing lepton. 
To ensure that the lepton charges are well measured, only events with |cos6'*| < 0.9 
are used. Furthermore, only unambiguously charged events are kept, i.e., the product 
of the charge of the two leptons must be —1. This requirement removes about 0.5% of 
the muon pairs and 6% of the tau pairs. The remaining charge misidentification level is 
estimated using simulated events to be 0.002% for muon pairs and 0.05% for tau pairs. 
The asymmetries are corrected for backgrounds and acceptance, using the appropriate 
Monte Carlo samples. This includes the effect of extrapolating from |cos^*| < 0.9 to 
Icos^l < 0.95. 



The measured and r+r^ asymmetries are shown in Table |T6|. For the Bhabha 

channel, rather than quoting an asymmetry, it is preferred to give the differential cross 
section with respect to 6*, because the reaction is dominated by t channel photon 



exchange. This is given in Table [17. 



A major contribution to the systematic uncertainties on the asymmetries is the 
background subtraction. The ditau channel is particularly sensitive to the Bhabha 
background since this is very peaked in the forward direction. For the tau pair asymmetry, 
this gives a correction of 1.3-3.8% for centre-of-mass energies of 130-183 GeV, respectively. 




Figure 10: Distribution of cos 9* in Bhabha events at 183 GeV with ^/s^/s > 
0.9. The data are represented by points, whilst the Monte Carlo expectation, 
normalized to the same integrated luminosity, is shown by the white histogram. 
The expected background is indicated by the shaded histogram. 




Figure 11: Distribution of cos 6'* in events at 183 GeV with ^Js'^/s > 

0.9. The data are represented by points, whilst the Monte Carlo expectation, 
normalized to the same integrated luminosity, is shown by the white histogram. 
The expected background is indicated by the shaded histogram. 



Table 15: Contributions to the systematic uncertainties on the ditau cross section 
measurements, for all energies and for both inclusive and exclusive processes. All 
quoted values are in percent. 





Description 


Eem (GeV) 


cut 




1 '^O 




1 fil 


1 72 


1 83 

± (JO 


0.1 


MC statistics 


n n 

u.y 


u.y 


1 1 

1.1 


1 1 
1.1 


l.U 




Energy scale 


1.0 


1.7 


1.0 


1.6 


1.5 




77 T^T~ 


0.5 


o.o 


0.7 


l.Z 


0.5 




77 


O.o 


0.9 


o.o 


0.9 


0.1 




qq 


0.4 


0.4 


0.1 


0.3 


0.3 




e+e~ 


n 7 


n 7 


u.o 


1 n 


u.o 




w+w- 








0.4 


0.7 




Luminosity 


1.0 


1.0 


0.7 


0.7 


0.5 


0.9 


MC statistics 


1.4 


1.3 


1.4 


1.2 


1.2 




Energy scale 


2.0 


2.0 


2.2 


2.1 


2.3 




77 — > T'^T~ 


0.5 


0.8 


0.6 


1.0 


0.4 




qq 


0.4 


0.6 


0.7 


0.3 


0.3 




e"''e~ 


1.3 


0.9 


1.0 


0.8 


0.9 






0.1 


0.1 


0.1 


0.1 


0.1 




w+w- 








0.2 


0.4 




Radiative background 


1.0 


1.4 


2.7 


0.5 


0.7 




Luminosity 


1.0 


1.0 


0.7 


0.7 


0.5 



The Monte Carlo statistical error on this correction enters as a part of the systematic 
uncertainty of the ditau asymmetry. Several other sources of systematic uncertainty are 
considered. The correction for event acceptance related to the extrapolation in polar angle 
introduces a significant uncertainty related to the finite Monte Carlo statistics. It also 
leads to a theoretical uncertainty due to ISR/FSR interference. This is assessed using 
ZFITTER. Systematic uncertainties from charge misidentification are negligible. 




cos(0*) 



Figure 12: Distribution of cos6'* in ditau events at 183 GeV with \/s'^/s > 
0.9. The data are represented by points, whilst the Monte Carlo expectation, 
normalized to the same integrated luminosity, is shown by the white histogram. 
The expected background is indicated by the shaded histogram. 



Table 16: Lepton forward-backward asymmetries with statistical and systematic 
uncertainties, calculated for ^/ s'/s >0.9 in the range |cos0| < 0.95. 





Lepton 


Afb 


SM prediction 


(GeV) 


Type 






130 




0.83 ± 0.08 ± 0.03 


0.70 






0.56 ± 0.12 ± 0.05 


0.70 


136 




0.63 ± 0.12 ± 0.03 


0.68 






0.65 ± 0.15 ± 0.04 


0.68 


161 




0.63 ± 0.11 ± 0.03 


0.61 






0.48 ± 0.14 ± 0.04 


0.61 


172 




0.72 ± 0.13 ± 0.04 


0.59 






0.44 ± 0.20 ± 0.04 


0.59 


183 




0.54 ± 0.06 ± 0.03 


0.58 






0.52 ± 0.08 ± 0.04 


0.58 



Table 17: Cross sections to produce electron pairs with \/s'/s > 0.9 and cos 6* 
in the quoted ranges. The quoted uncertainties include statistical and systematic 
components. 



Eem (GeV) 


cos ^min; 




a (pb) 


SM prediction 


130 


-0.9, 


-0.7 


0.19 ± 0.34 


0.37 




-0.7, 


-0.5 


1.41 ± 0.35 


0.55 




-0.5, 


-0.3 


1.36 ± 0.45 


1.09 




-0.3, 


-0.1 


1.23 ± 0.48 


1.19 




-0.1, 


0.1 


2.60 ± 0.69 


2.45 




0.1, 


0.3 


3.78 ± 0.83 


3.82 




0.3, 


0.5 


8.88 ± 1.18 


7.36 




0.5, 


0.7 


21.63 ± 2.12 


22.20 




0.7, 


0.9 


149.61 ± 6.22 


148.00 


136 


-0.9, 


-0.7 


0.73 ± 0.20 


0.22 




-0.7, 


-0.5 


1.16 ± 0.36 


0.62 




-0.5, 


-0.3 


0.54 ± 0.35 


0.49 




-0.3, 


-0.1 


0.52 ± 0.41 


0.89 




-0.1, 


0.1 


1.46 ± 0.62 


2.09 




0.1, 


0.3 


2.09 ± 0.74 


2.96 




0.3, 


0.5 


6.68 ± 1.08 


6.13 




0.5, 


0.7 


16.58 ± 1.97 


20.50 




0.7, 


0.9 


132.55 ± 5.85 


133.00 


161 


-0.9, 


-0.7 


0.46 ± 0.21 


0.37 




-0.7, 


-0.5 


0.88 ± 0.21 


0.44 




-0.5, 


-0.3 


0.55 ± 0.28 


0.79 




-0.3, 


-0.1 


0.39 ± 0.26 


0.62 




-0.1, 


0.1 


1.24 ± 0.40 


1.43 




0.1, 


0.3 


2.37 ± 0.47 


2.07 




0.3, 


0.5 


5.35 ± 0.73 


4.95 




0.5, 


0.7 


14.38 ± 1.27 


14.10 




0.7, 


0.9 


93.76 ± 3.76 


94.20 


172 


-0.9, 


-0.7 


0.32 ± 0.19 


0.28 




-0.7, 


-0.5 


0.88 ± 0.19 


0.34 




-0.5, 


-0.3 


0.66 ± 0.24 


0.58 




-0.3, 


-0.1 


0.61 ± 0.23 


0.44 




-0.1, 


0.1 


0.95 ± 0.36 


1.23 




0.1, 


0.3 


1.80 ± 0.47 


1.93 




0.3, 


0.5 


4.92 ± 0.71 


4.24 




0.5, 


0.7 


13.07 ± 1.20 


12.40 




0.7, 


0.9 


84.61 ± 3.51 


81.10 


183 


-0.9, 


-0.7 


0.24 ± 0.07 


0.21 




-0.7, 


-0.5 


0.29 ± 0.07 


0.25 




-0.5, 


-0.3 


0.46 ± 0.10 


0.51 




-0.3, 


-0.1 


0.71 ± 0.12 


0.64 




-0.1, 


0.1 


0.83 ± 0.14 


0.90 




0.1, 


0.3 


1.42 ± 0.20 


1.83 




0.3, 


0.5 


3.90 ± 0.29 


3.66 




0.5, 


0.7 


12.47 ± 0.56 


11.10 




0.7, 


0.9 


71.90 ± 1.86 


71.80 



7 Interpretation in Terms of New Physics 



7.1 Comparison with Standard Model Predictions 

The measured cross sections and asymmetries from Tables ^ and |16| are plotted as 



a function of centre-of-mass energy in Figs. ^ and [IJ, respectively. The results are 
compared with SM predictions based on BHWIDE |^ for electron pair production and 
ZFITTER 1 21 1 for all other processes. The ZFITTER predictions [| are computed from 
the input values mz = 91.1867 GeV/c^, = 174.1 GeV/c^, mn = 127.0 GeV/c^, 
aem(^z) = 1/128.896 and as{Mz) = 0.120. The exclusive cross sections and asymmetries 
are given in the restricted angular range for the outgoing fermion direction, |cos^| < 0.95. 
The inclusive cross sections correspond to the full angular range. 
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Figure 13: Measured cross sections for fermion pair production. The curves 
indicate the predictions obtained from BHWIDE for the Bhabha process and 
from ZFITTER for the other channels. (Some of the points are shifted slightly 
along the horizontal axis to prevent them overlapping). 



Because of the poor knowledge of the contribution of ISR/FSR interference, the 
ZFITTER predictions are assigned a systematic uncertainty equal to the difference in 
the predictions with and without interference. This amounts to 1.5% for the hadronic 
cross section and 2% for the t~^t~ , bb and cc cross sections. These uncertainties 

would almost double if one attempted to extrapolate the results to the full angular range. 
The Bhabha cross section is assumed to be uncertain by 3%, based on a study of different 

1 Default ZFITTER flags are used, except for BOXD = 1, CONV = 1, INTF = 1 and INCL = 0. The 
flag FINR = is used for hadronic events and FINR = 1 for dilepton events. 
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Figure 14: Measured asymmetries for muon and tau pair production. The curves 
indicate the predictions obtained from ZFITTER. (Some of the points are shifted 
slightly along the horizontal axis to prevent them overlapping). 



event generators. These uncertainties are taken into account in the calculation of limits 
on physics beyond the SM given in the following sections. 

The measured cross sections and asymmetries are consistent with SM predictions. 
Similar results have been published by the OPAL Collaboration The L3 



Collaboration has published results covering centre-of-mass energies up to 172 GeV |23] 



7.2 Limits on Four-Fermion Contact Interactions 

Comparing the measured exclusive difermion cross sections and angular distributions 
with SM predictions allows one to place limits upon many possible extensions to the 
SM. One convenient parametrization of such effects is given by the addition of four- 



fermion contact interactions to the known SM processes. Such contact interactions 
are characterized by a scale A, interpreted as the mass of a new heavy particle exchanged 
between the incoming and outgoing fermion pairs, and a coupling g giving the strength 
of the interaction. Contact interactions are, for example, expected to occur if fermions 
are composite. 



Following the notation of Ref. |2^, the effective Lagrangian for the four- fermion 



contact interaction in the process e^e — > ff is given by 



l^i -I- 0)1\ ij=L,R 



with 5 = 1 if f = e, or otherwise. The fields ei^pt (fi,/?) are the left- and right-handed 



chirality projections of electron (fermion) spinors. The coefficients rjij, which take a 
value between —1 and +1, indicate the relative contribution of the different chirality 
combinations to the Lagrangian. The sign of ?7sign determines whether the contact 
interaction interferes constructively or destructively with the SM amplitude. Several 
different models are considered in this analysis, corresponding to the choices of the resign 



and rjij given in Table 18 



Table 18: Four-fermion interaction models. 



Model 


'7sign 


Vll 


Vrr 


'7LR 


'7RL 


LL± 


±1 


1 











RR± 


±1 





1 








VV± 


±1 


1 


1 


1 


1 


AA± 


±1 


1 


1 


-1 


-1 


LR± 


±1 








1 





RL± 


±1 











1 


LL+RR± 


±1 


1 


1 








LR+RL± 


±1 








1 


1 



In the presence of contact interactions the differential cross section for e"'"e~ —>■ ff as 
a function of the polar angle 6 of the outgoing fermion with respect to the e~ beam line 
can be written as 



d(J „ , , 

'^SM(s,t) 



dcosO 



(10) 



with s and t being the Mandelstam variables and e = g'^r]sign/ {4:7rA'^). The SM cross 
section Fsm is computed as described in Section fTII . The contributions to the cross section 
from the SM - contact interaction interference and from the pure contact interaction are 
denoted by F^°'^'^ and F^^^^, respectively. They are calculated in the improved Born 



approximation. The Born level formulae can be found in Ref. ||2^ and these are corrected 



for ISR according to Ref. |^6[. Because no higher order calculations are available for the 
contact interactions, the ratios of these with the improved Born predictions for the SM 
cross sections are taken, to allow for a partial cancellation of higher order effects. 

The predictions of Equation |10| are fitted to the data using a binned maximum 
likelihood method. For contact interactions affecting the dilepton channels, the likelihood 
function £ is defined by 

r /^^^corr. T\ I I /^/^uncorr. -i \ I I -7-) / Ardata lArPred/ \ , corr A„corr 1 ^ uncorr A „ uncorr \ 

L = G{a ,l)[[G{ai A) [[V [N^^ , [iVf^ (e) + a Ara-^ + IXn,-^^ Jj. 

i k 

(11) 

The indices i and k run over the centre-of-mass energy points and angular bins in cos^*, 
respectively. The function V gives the Poisson probability to observe Nf^'''^ events in 
the data if iV^^'^*''^ are expected. The systematic uncertainties on the expected number of 
events which are (un) correlated between the centre-of-mass energy points are represented 
by (An'^"™") An™", respectively. These uncertainties are taken into account using the 
parameters a^°^^ and a""™'''', which are constrained using Gaussian distributions G with 
zero mean and unit standard deviation. These parameters are fitted together with the 



parameter e. The fit range in the angular distribution is chosen to be |cos6'*| < 0.9 for 
the Bhabha channel and |cos^*| < 0.95 for the muon and tau pair channels. 

For contact interactions affecting hadronic events, the sum over angular bins is 
dropped, and instead two additional terms are added to the likelihood function to take 
into account the constraints from the measurement of the jet charge asymmetries, given 
by Equation | and Table |ll]. The contact interaction can be assumed to couple to all 
quark flavours with equal strength. In this case, the jet charge asymmetry measurements 
improve the limits for some models by up to 70%, primarily because of their sensitivity 
to the relative cross sections of up- and down-type quarks. Alternatively, one can assume 
that contact interactions only affect the bb final state. Such interactions are strongly 



constrained using the measurement of -Rb from Section |5]^ and the jet charge asymmetry 
in b-enriched events of Section |5.4| . 

Because of the quadratic dependence of the theoretical cross sections upon e, the 
likelihood function can have two maxima. The 68% confidence level limits (e"*" and e~) 
on e are therefore estimated as follows: 



/6 /"oo /"oo 

C{e')dt' = / £(e')c/e' = 0.16 / C{t')de' 



(12) 



where for each value of e' the parameters ac and at are chosen which maximize the 
likelihood. The results for contact interaction affecting leptonic final states are listed 
in Table ITPI. Table EUI gives the results obtained for contact interactions affecting both 



hadronic events and all difermion events. 

Although all the physics content is described by the well-defined parameter e, it is 
conventional to extract limits on the energy scale A, assuming = An. The 95% 
confidence level limits €95 are computed according to 

Z^"' C{e')de' = 0.95 H C{e')de' , /° C{e')de' = 0.95 T C{e')de' , (13) 

Jo Jo JSa^ J-OQ 



which are then used to obtain 



A^ 



^95 I 



(14) 



Limits on the energy scale are listed in Tables and One can drop the assumption 
g"^ = 4:71, in which case these results become limits on y/AirA/g. 

These results are competitive with previous analyses of contact interactions already 
performed at LEP |22| , at the Tevatron p8| and at HERA ||29|. However, models of 



e^e~uii and e^e~dd contact interactions which violate parity (LL, RR, LR and RL) are 
already severely constrained by atomic physics parity violation experiments, which quote 
limits of the order of 15 TeV ||3^. The LEP limits for the fully leptonic couplings or 
those involving b quarks are of particular interest since they are inaccessible at pp or ep 
colliders. 



7.3 Limits on R-parity Violating Sneutrinos 

Supersymmetric theories with R-parity violation have terms in the Lagrangian of the form 
XijkLiLjEk, where L denotes a lepton doublet superfield and E denotes a lepton singlet 
superfield. The parameter A is a Yukawa coupling, and i, j, k = 1, 2, 3 are generation 
indices. The Xijk, which for the purposes of this analysis are assumed to be real, are 



Table 19: Limits on contact interactions coupling to dilepton final states. The 
68% confidence level range is given for e whilst the 95% confidence level limits are 
given for A. The results presented for l+l" assume lepton universality. 



iVlOQcl 


F >^ J l-"-®^ ) 


i V ±e V j 


J V xe V j 


+ — + — 

e^e — s> e^e 








T T 

LL 


[— U.UO 1, +U.UziJ 


6.1 


0.0 


T-> 


[— U.Ud7, +U.U2/J 


6.1 


3.4 


V V 


[—0.017, +U.UUc)J 


6.4 


8.0 


A A 

AA 


[— (J.U18, +0.U19J 


4.2 


5.5 


T T ) 


[—0.042, +0.015J 


A A 

4.0 


A 

4.2 


T T 1 T) T) 

LL+KK 


r A A O O 1 A A A A1 

[— O.Odo, +0.009J 


4.2 


r A 

5.0 


LK+RL 


r A AOO 1 A AA/^1 

[—0.022, +0.006J 


5.5 


6.5 


+ — + — 
e^e — > ijl^h 








T T 

LL 


r nni/i inn/inl 
[—0.014, +0.040J 


A '7 

4. i 


4.0 


RR 


r AA1/^ iAA/1 Ol 

[—0.016, +0.043J 


A A 

4.4 


3.8 


V V 


r A AAr 1 A A1 

[—0.005, +O.OI0J 


1.1 


6.3 


A A 

AA 


r A AAA 1 A A1 rl 

[—0.009, +0.015J 


6.8 


6.Z 


LR 


r A O^A 1 A AOrl 

[—0.270, +0.025J 


1.8 




3.8 


T T 1 T) T) 

LL+RR 


r A AA'7 1 A AOOl 

[—0.007, +0.02zJ 


6.6 


5.4 


LR+RL 


r A C\nc\ 1 A A"! A1 

[-0.260, +0.019J 


1 A 

1.9 


5.1 


-I- — -t- — 
e^e — >• T^T 








T T 

LL 


[ — U.Uoy, -hU.Uozj 


Q 7 
0. / 


o.y 


RR 


[-0.046, +0.034] 


3.4 


3.7 


VV 


[-0.012, +0.016] 


6.2 


5.9 


AA 


[-0.022, +0.013] 


5.2 


5.6 


LR 


[-0.275, +0.033] 


1.8 


3.3 


LL+RR 


[-0.020, +0.018] 


5.2 


5.2 


LR+RL 


[-0.265, +0.025] 


1.8 


4.3 


e+e- ^ 1+1- 








LL 


[-0.014, +0.020] 


5.5 


5.3 


RR 


[-0.016, +0.021] 


5.3 


5.1 


W 


[-0.005, +0.006] 


9.5 


9.3 


AA 


[-0.007, +0.010] 


8.0 


7.5 


LR 


[-0.023, +0.019] 


4.8 


5.0 


LL+RR 


[-0.008, +0.010] 


7.7 


7.3 


LR+RL 


[-0.011, +0.009] 


7.1 


7.2 



Table 20: Limits on contact interactions coupling to hadronic or to all difermion 

final states. The 68% confidence level range is given for e whilst the 95% confidence 
level limits are given for A. The results presented for ff assume that the contact 
interaction couples to all the outgoing fermion types equally. 
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non- vanishing only for i < j. These terms allow for single production of sleptons at e^e~ 
collider experiments. 

At LEP, dilepton production cross sections could then differ from their SM 
expectations as a result of the exchange of R-parity violating sneutrinos in the s or 



t channels |^T]]. Table |21] shows the most interesting possibilities. Those involving 
s channel sneutrino exchange lead to a resonance. For the results presented here, this 
resonance is assumed to have a width of 1 GeV/c^, which can occur if the sneutrino 
also has R-parity conserving decay modes For a sneutrino only having the R-parity 
violating decay mode into lepton pairs, the width would be much less than this, leading 
to slightly better limits. 



Table 21: For each dilepton channel, the table shows the coupling amplitude, 
the sneutrino type exchanged, and an indication of whether the exchange occurs 
in the s or t channel. 
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Direct searches for R-parity violating sneutrinos at LEP have led to lower limits on 
their masses of 72 GeV/c^ for Uf, and 49 GeV/c^ for z/^ and Ur Indirect limits based 
upon lepton universality and leptonic tau decays imply that the couplings in Table pi] 
must approximately satisfy A < 0.1(m;-^/200) GeV/c^, where mj^ is the mass of the 
appropriate right-handed charged slepton These limits can be improved using the 



dilepton cross section and asymmetry data presented in this paper. 

Limits on the couplings are obtained by comparing the measured dilepton differential 
cross sections with respect to the polar angle with the theoretical cross sections in reference 
|]3T| . The likelihood function used in the fit and the corrections for ISR, etc. follow the 
procedure of Section |7.2|. The fit is performed in terms of the parameter A^. Since the 



likelihood function can have two minima, limits are again determined by integrating the 
likelihood function with respect to A^. A one-sided limit is used when A^ is positive 
definite, which occurs when A^ = A^2i -^1315 ^'^^ when A^ = A121A233 or Ai3iA232- 



Figures |T5|, and |1^ show the results for those processes involving sneutrino exchange 



in the s channel. Similar results have been obtained by the OPAL Collaboration P2[ and 



for Ur by the L3 Collaboration |33 



Limits on A121 and A131 from t channel exchange of in muon and tau pair production, 
respectively, give much weaker limits. These rise from |Aiji| < 0.5 at z/g = 100 GeV/c^ to 
< 0.9 at Ue = 300 GeV/c^. 

7.4 Limits on Leptoquarks and R-Parity Violating Squarks 

At LEP, the t channel exchange of a leptoquark can modify the qq cross section and 



jet charge asymmetry, as described by the Born level equations given in Ref. |3^. A 
comparison of the measurements with these equations allows upper limits to be set on the 
leptoquark's couplings g as a function of its mass Mlq, using the same fit technique and 
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Figure 15: 95% confidence level upper limits, obtained from the Bhabha cross 
sections, on |Ai2i| versus the assumed mass and on |Ai3i| versus the assumed 
z>T- mass. 
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Figure 16: 95% confidence level upper limits, obtained from the /x^/x cross 
sections, on vT^i3i^^232l versus the assumed mass. 




Figure 17: 95% confidence level upper limits, obtained from the ditau cross 
sections, on ^|Ai2iA233| versus the assumed t'^ mass. 



corrections for ISR, etc. as employed in Section 7^. Although the leptoquark t channel 
exchange alters the angular distribution of the outgoing qq system, this has negligible 
effect on the qq selection efficiency. Limits are obtained for each possible leptoquark 
species. The allowed species can be classified according to their spin, weak isospin I 
and hypercharge. Scalar and vector leptoquarks are denoted by symbols S/ and V/ 
respectively, and isomultiplets with different hypercharges are distinguished by a tilde. 
An indication "(L)" or "(R)" after the name indicates if the leptoquark couples to left- or 
right-handed leptons. Where both chirality couplings are possible, limits are set for the 
two cases independently, assuming that left- and right-handed couplings are not present 
simultaneously. It is also assumed that leptoquarks within a given isomultiplet are mass 
degenerate . 

The Si (L) and So (L) leptoquarks are equivalent to up-type anti-squarks and down- 
type squarks, respectively, in supersymmetric theories with an R-parity breaking term 
X[jiJ-iiQjT)k {j,k = 1,2,3). Limits in terms of the leptoquark coupling are then exactly 
equivalent to limits in terms of A'^^,. 

Table ^ gives for each leptoquark type, the 95% confidence level lower limits on its 
mass Mlq, assuming that it has a coupling strength equal to the electromagnetic coupling 
g = e. The limits are given separately, assuming that (i) the leptoquark couples to only 
]^st only 2'^'^ generation quarks, or (ii) to only 3'''^ generation quarks. The former limits 
are derived using the measured hadronic cross section and jet charge asymmetries, whilst 
the latter uses the bb cross section and jet charge asymmetries. For g ^ e, the mass limit 
scales approximately in proportion to the coupling if it exceeds about 200 GeV/c^. (This 
is the contact term limit.) 

Similar limits have been obtained by OPAL and L3 [E^, 12^. Limits from the 



Tevatron |^ |3^ depend upon the assumed branching ratio of the charged leptonic decay 



Table 22: 95% confidence level lower limits on the leptoquark mass for each 
species. Limits are given separately according to the quark generation to which 
the leptoquark is assumed to couple. A dash indicates that no limit can be set 
and "N.A." denotes leptoquarks coupling only to top quarks and hence not visible 
at LEP. 



Quark 
Generation 


Limit on scalar leptoquark mass (GeV/c^) 


So (L) So (R) So (R) Si (L) Si (R) Si (L) Si (L) 


pt 

3rd 


200 — 70 240 — 20 — 
N.A. N.A. 180 450 — N.A. 50 



Quark 
Generation 


Limit on vector leptoquark mass (GeV/c^) 


Vi (R) Vi (L) Vi (L) Vo (L) Vo (R) Vo (R) Vi (L) 


3rd 


150 130 90 340 120 280 470 
260 160 N.A. 400 140 N.A. 400 



mode. However, if this is 100%, then the Tevatron excludes leptoquark masses below 
~ 225 GeV/c^. Other experiments can place limits on leptoquarks which couple to first 
generation quarks. In particular, low energy data such as atomic parity violation and rare 
decays give very stringent limits, usually as a function of the ratio ^^/Mlq. If g = e, they 
imply a lower limit on the leptoquark mass in the range 430-1500 GeV/c^ |3^, depending 
on the leptoquark species. Preliminary results from HERA |3^ exclude scalar leptoquarks 
with masses below ~ 250 GeV/c^ if g = e. 



7.5 Limits on Extra Z Bosons 

To unify the strong and electroweak interactions. Grand Unification Theories (GUT) 
extend the SM gauge group to a group of higher rank, predicting therefore the presence of 
at least one extra neutral gauge boson Z'. The theories which are considered in this section 
are the Eq and the Left-Right (LR) models . In the Eq model, the unification group 
Eq can break into the SM SUciS) ® SUl{'2^) ® ^y(l) in different ways. Each symmetry 
breaking pattern leads to the presence of at least one extra U(l) symmetry and therefore 
one extra gauge boson. This is characterized by a parameter 6e(, (— 7r/2 < 6'£;g< 7r/2) 
which entirely defines its couplings to conventional fermions. Contributions from exotic 
particles predicted by E^ or supersymmetric particles are ignored in this analysis. Four 
models derived from Eq are studied here, the Eq{x)i Eq{iP), Eq^t]) and Eq{I) defined by 



(^Ee= 0, 7r/2, — arctany5/3 and arctany3/5, respectively In the LR model, the 

SM group is extended to SUc{3) (g) SUl{2) O SUr{2) O Ub-l{1), where B and L are 
the baryon and lepton number, respectively. The SM f/y(l) symmetry is recovered by 
a linear combination of the generator of Ub~l{^) and the third component of SUr{2). 
This model, which can arise from symmetry breaking of the GUT group SO(IO), leads 
to the presence of one extra neutral gauge boson. Contributions from the extra of 



SU{2)[i are neglected in this analysis. Couplings to conventional fermions depend only 



on one parameter aiR, where y 2/3 < tti/j < y (1 — 2 sin 6'iy)/sin 6w, which is function 
of the couphng constants gL,R of SUl,r{2) and sin^6'iy. Only the Left-Right Symmetric 
Model (LRS) will be studied here, which has ql = Qr, implying that if sin^ 9]^ = 0.23 
then aLR = 1.53 . 

Outside the context of GUT, a limit is derived on the mass of the non-gauge-invariant 
sequential SM (SSM) Z' boson, having the same couplings as the SM Z, but with a higher 
mass. Limits are also placed upon the axial and vector couplings of an arbitrary Z' as a 
function of its mass. 

In all the models mentioned above, the symmetry eigenstate Z'° of the extra U(l) or 
SUr{2) can mix with the symmetry eigenstate Z° of SUl(2) f/y(l) with a mixing angle 
^mix- In such a case, the Z resonance observed at LEP I must be identified as one of the 
mass eigenstates of the Z'^-Z^ system while the second mass eigenstate Z' of mass mz' is 
a free parameter ^ | . 

To obtain the 95% confidence level exclusion limits on the various free parameters, 
least squares fits are performed using the set of ALEPH measurements given in Table 23, 
taking into account the correlations between them. The LEPl measurements are taken 
from Ref. |E2|, whilst the LEP2 measurements are presented in Tables 0, 1^ and |16. 



Table 23: Set of observables used in the Z' analyses and the corresponding SM 
values per degree of freedom. 





Observables 


xIm/ndf. 


LEPl 


" 5 ^FB 1 ^ 1 — h^i' 


130.9/120 


LEP2 
130-183 GeV 


/1 1+1- _qq D 

1 = fi,T, ^Js'/s > 0.9, Icos^l < 0.95 


19.4/29 



Theoretical predictions for difermion cross sections and asymmetries are obtained from 
the program ZEFIT 5.0 [^, which is an extension of ZFITTER 5.0 including the 
one extra neutral gauge boson of the Eq, LR or SSM models. Theoretical uncertainties 
on the SM predictions are taken into account as described in Section |7]T|. For all models, 
the minimum is found to occur when no Z' boson is present. 

For the five models Eeix), E^i^), Ee{T]), Eq{I) and LRS, Fig. |TB| shows the 95% 
confidence level limits obtained in the plane of Z' mass versus mixing angle 6'mix. Both 
parameters are treated as independent, so these limits correspond to a increase of 
5.99 . The LEPl data mainly constrain the mixing angle, whilst the LEP2 data mainly 
constrain the Z' mass at small mixing angles. 

Alternatively, assuming ^mix = 0, lower limits on the Z' mass can be obtained using 
a one-sided, one-parameter fit (Ax^ = 2.71). The resulting limits are given in Table |2^ , 
where they are compared with those from direct Z' searches performed by the CDF 
Collaboration [^. This table also gives the mass limit for the SSM Z', which is superior 
to the limit from CDF. 

Limits can also be placed on the vector and axial couplings of an arbitrary Z', as 
a function of its mass. To simplify, such limits will only be given here for the leptonic 
couplings (assuming lepton universality) and also only for the case 6'jnix = 0. Limits are 
placed on the two couplings simultaneously (Ax^ = 5.99). The excluded region is found 



Table 24: Comparison of 95% confidence level lower limits on niz' (GeV/c^) from 
one parameter electroweak fits (ALEPH) and direct searches (CDF) for ^mix = 0. 
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Figure 18: 95% confidence level limits in the niz' vs. ^mix plane for the Eq(x)j 
Eq{iI;), E^iri), Eq{I) and LRS models. The shaded regions are excluded. 



to be approximately rectangular in shape and its size is given as a function of Z' mass in 
Table 



Table 25: 95% confidence level limits on the axial g'^ and vector g'^ couplings of 
a Z' boson of mass mz' to a lepton pair. 



mz' (GeV/c2) 




9'a 






9'. 




300 
600 
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< 1.11 



8 Conclusions 

Measurements of the hadronic and leptonic cross sections and asymmetries at y/s = 
130-183 GeV have been presented. The ratios of the bb to qq production cross sections 
at y/s = 130-183 GeV and of the cc to qq production cross sections at y/s = 183 GeV have 
been shown, as well as jet charge asymmetries. The results agree with the predictions of 
the Standard Model and allow limits to be placed on four-fermion contact interactions, 
R-parity violating sneutrinos and squarks, leptoquarks and Z' bosons. The limits on 
the energy scale A of e"'"e~ff contact interactions are typically in the range from 2- 
10 TeV. Those for e"'"e~l"'"l~ and e"'"e~bb interactions are of particular interest, since 
they are inaccessible at colliders using proton beams. The new ALEPH limits on R- 
parity violating sneutrinos reach masses of a few hundred GeV/c^ for large values of their 
Yukawa couplings. 
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